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We analyze dilepton emission from hot and dense matter using a hybrid approach based on the 
Ultrarelativistic Quantum Molecular Dynamics (UrQMD) transport model with an intermediate 
hydrodynamic stage for the description of heavy-ion collisions at relativistic energies. During the 
hydrodynamic stage, the production of lepton pairs is described by radiation rates for a strongly 
interacting medium in thermal equilibrium. In the low mass region, hadronic thermal emission 
is evaluated assuming vector meson dominance including in-medium modifications of the p meson 
spectral function through scattering from nucleons and pions in the heat bath. In the intermedi- 
ate mass region, the hadronic rate is essentially determined by multi-pion annihilation processes. 
Emission from quark-antiquark annihilation in the quark gluon plasma is taken into account as well. 
When the system is sufficiently dilute, the hydrodynamic description breaks down and a transition 
to a final cascade stage is performed. In this stage dimuon emission is evaluated as commonly done 
in transport models. Focusing on the enhancement with respect to the contribution from long-lived 
hadron decays after freezout observed at the SPS in the low mass region of the dilepton spectra, 
the relative importance of the different thermal contributions and of the two dynamical stages is 
investigated. We find that three separated regions can be identified in the invariant mass spec- 
tra. Whereas the very low and the intermediate mass regions mostly receive contribution from the 
thermal dilepton emission, the region around the vector meson peak is dominated by the cascade 
emission. Above the p-peak region the spectrum is driven by QGP radiation. Analysis of the dimuon 
transverse mass spectra reveals that the thermal hadronic emission shows an evident mass ordering 
not present in the emission from the QGP. A comparison of our calculation to recent acceptance 
corrected NA60 data on invariant as well as transverse mass spectra is performed. 
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I. INTRODUCTION 

Dileptons represent a penetrating probe of the hot and 
dense nuclear matter created in heavy ion collisions since, 
once produced, they essentially do not interact with the 
surrounding hadronic matter. The analysis of the electro- 
magnetic response of the dense and hot medium is tightly 
connected to the investigation of the in-medium modifi- 
cation of the vector meson properties. Vector mesons are 
ideally suited for this exploration, because they can di- 
rectly decay into a lepton-antilepton pair. One therefore 
aims to infer information on the modifications induced 
by the medium on specific properties of the vector me- 
son, such as its mass and/or its width, from the invariant 
mass dilepton spectra. 

The NA60 experiment at CERN SPS has recently mea- 
sured the invariant mass spectra of low mass dimuons 
3] . As in previous ultra-relativistic heavy ion collision ex- 
periments [1, @] , at low invariant masses an enhancement 
of dilepton pairs in heavy systems as compared to the 
expected contribution from decays of long-lived hadrons 
(mostly 77, 77', uj and 4> mesons), commonly referred to as 
"hadron cocktail" , was observed. The high data quality 
of NA60 allowed to subtract the cocktail of the decay 
sources from the total data and to isolate the excess of 
pairs. One of the first results was a strong evidence for 



the presence of broadening of the spectral function of the 
p meson 

From the first communication by the NA60 collabo- 
ration m, the remaining excess of pairs (often referred 
to as "the excess") has been object of various theoreti- 
cal investigations. Phenomenological models have been 
refined and extended to include the various sources of 
dilepton production, and an almost comprehensive 
interpretation of the invariant mass spectra of the excess 
has been reached. Calculations performed by convolu- 
tion of dilepton production rates, using the assumption 
of local equilibrium, over the space-time evolution of the 
medium modelled according to simple expanding fireball 
approaches have been quite successful in explaining the 
dilepton excess observed at SPS 7, 9, 11, 12]. 

This use of basic approaches for the modelling of the 
heavy-ion collision dynamics focused the investigations 
on the various scenarios for the in-medium modifica- 
tion of the vector meson properties, and an evidence for 
a broadening of the p meson spectral function was in- 
ferred from the comparison to the measurements. How- 
ever, in these models the time development of the reac- 
tion is not described by fully dynamical simulations but 
rather parametrized in terms of the estimated time which 
the system spends in the various phases. Though such 
parametrizations can be properly tuned and constrained. 
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a more complete treatment of the heavy-ion dynamics in 
terms of transport and/or fuUy (3+1) dimensional hydro- 
dynamical models should be in order, since a good knowl- 
edge of the temperature, baryon-density and flow evolu- 
tion is necessary to convert dilepton rates into space-time 
integrated spectra. 

Earlier, a (2-f 1) dimensional hydrodynamical descrip- 
tion of dilepton production has been performed by Huo- 
viven et al. jl3l - [l5j and compared to data from the 
first generation of dilepton experiments performed in the 
nineties. However, already at that time [15|| it was shown 
that the evolution can have strong impact on the dilepton 
yields. To our knowledge, the model was never applied to 
the high resolution NA60 data. Application of hydrody- 
namics to analyse the dimuon excess observed by NA60 
were performed by Dusling et al. [1, [13] using a boost 
invariant (1+1) hydrodynamical model. If on one side 
this represents a step further towards the inclusion of a 
more realistic dynamics, on the other side boost invari- 
ant hydrodynamics is expected to be an approximation 
more valid at RHIC energies than at SPS. In Ref. [l6l |. 
e.g., it was shown that the photon spectra measured by 
the WA98 Collaboration at the CERN SPS ^ could 
be accounted for in a boost-invariant scenario only by 
choosing a short initial time (ro=0.3 fm) which is not in 
accordance with the expected longitudinal geometry. In 
this respect, it is surely desirable to relax this approxi- 
mation and use a full (3+1) hydrodynamic expansion. 

Moreover, the contribution to the excess of non- 
thermal dilepton radiation from short living mesons has 
only very recently received proper attention by the theo- 
retical groups mostly active in the field of thermal dilep- 
tons. For example, in Ref. [8j] the contribution from 
freezeout p mesons was neglected and added only in a 
later work [l3|, whereas in Ref. [l^] an explicit treat- 
ment of p decays at thermal freezeout together with the 
additional inclusion of a non-thermal component of pri- 
mordially produced p mesons which escape the fireball 
were required in order to account for some discrepancy 
with thermal emission revealed by analysis of measured 
transverse pair momenta (pr) spectra at pt>1 GeV. 

In this paper, we want to directly address these aspects 
presenting a consistent calculation of the dilepton pro- 
duction at SPS energy within a model which attempts 
to take into account both the complexity of the dilep- 
ton rate in hot a dense medium as well as the complex- 
ity of the pre-, post-, and equilibrium heavy-ion dynam- 
ics. The latter is modelled with an integrated Boltz- 
mann+hydrodynamics hybrid approach based on the Ul- 
trarelativistic Quantum Molecular Dynamics (UrQMD) 
transport model with an intermediate (3+1) dimensional 
ideal hydrodynamic stage. For earlier investigations 
within hybrid approaches see e.g. [l8l - |20j . The hybrid 
approach used here has been successfully applied to many 
bulk observables and photon spectra [21| at SPS energies. 
For the present investigation, during the local equilib- 
rium phase, the radiation rate of the strongly interacting 
medium is standardly modelled using the vector meson 



dominance model and related to the spectral properties 
of the light vector mesons, with the p meson having the 
dominant role. In-medium modifications of the p-meson 
spectral function due to scattering from hadrons in the 
heat bath are properly included in the model. Two addi- 
tional sources of thermal radiation, namely emission from 
multi-pion annihilation processes and from a thermalized 
partonic phase are included as well. 

The paper is organized as follows. In Sec. |TT1 we briefly 
discuss the hybrid model and present the emission rates 
of the various sources taken into account. In Sec. IIIII 
we present results for the dilepton excess invariant mass 
spectra and transverse mass spectra. A comparison to 
the NA60 data is performed. Section IIVI is dedicated to 
a discussion of the the various components that enter the 
contribution of the cascade stage of the hybrid approach. 
Finally, a summary and conclusions are given in Sec. [V] 



II. THE MODEL 
A. The hybrid approach 

To simulate the dynamics of the In+In collisions we 
employ a transport approach with an embedded three- 
dimensional ideal relativistic one fluid evolution for the 
hot and dense stage of the reaction based on the UrQMD 
model. The present hybrid approach has been exten- 
sively described in Ref. [1^ and first a ppli cations to e.m. 
probes have been recently performed [2ll [23l - l26| . Here, 
we limit ourselves to briefly describe its main features 
and refer the reader to Ref. for details. 

UrQMD [27l - l29j is a hadronic transport approach 
which simulates multiple interactions of ingoing and 
newly produced particles, the excitation and fragmen- 
tation of color strings and the formation and decay of 
hadronic resonances. The coupling between the UrQMD 
initial state and the hydrodynamical evolution proceeds 
when the two Lorentz-contracted nuclei have passed 
through each other. Here, the spectators continue to 
propagate in the cascade and all other hadrons are 
mapped to the hydrodynamic grid. Event-by-event fluc- 
tuations are directly taken into account via initial con- 
ditions generated by the primary collisions and string 
fragmentations in the microscopic UrQMD model. This 
leads to non-trivial velocity and energy density distribu- 
tions for the hydrodynamical initial conditions in each 
single event [30, |31|. Subsequently, a full (3+1) dimen- 
sional ideal hydrodynamic evolution is performed using 
the SHASTA algorithm jH,!!!. The hydrodynamic evo- 
lution is gradually merged into the hadronic cascade: 
to mimic an iso-eigentime hypersurface, full transverse 
slices, of thickness Az — 0.2fm, are transformed to par- 
ticles whenever in all cells of each individual slice the en- 
ergy density drops below five times the ground state en- 
ergy density. The employment of such gradual transition 
allows to obtain a rapidity independent transition tem- 
perature without artificial time dilatation effects [13, [IH] 



3 



and has been explored in detail in various recent works 
[3^ [s^-IHl devoted to SPS conditions. When merging, 
the hydrodynamic fields are transformed to particle de- 
grees of freedom via the Cooper-Frye equation [s^ on the 
hypersurface cr^ 



= / f{x,p)p''daf. 



(1) 



where f{x,p) are the boosted Fermi or Bose distributions 
corresponding to the respective particle species. For the 
present analysis, the latter was modified in order to ac- 
count for the spectral shape of the p meson via the sub- 
stitution: 



d^p d^p 



Po 



Po 



dM^S+ (M^ - ml) d*p2Ap{M) , (2) 



with Ap{M) = — ilmDp(M) being the p meson spec- 
tral function and Dp{M) being the p meson propagator 
in vacuum. The created particles then proceed in their 
evolution in the hadronic cascade where rescatterings and 
final decays occur until all interactions cease and the sys- 
tem decouples. 

An input for the hydrodynamical calculation is the 
equation of state (EoS). In this work we employ an equa- 
tion of state in line with lattice data that follows from 
coupling the Polyakov loop to a chiral hadronic flavor- 
SU(3) model [13] • The hadronic part is an extension of a 
non-linear representation of a sigma-omega model includ- 
ing the lowest-lying multiplets of baryons and mesons (for 
the derivation and a detailed discussion of the hadronic 
part of the model Lagrangian see [4]1443| ). In spirit simi- 
lar to the PNJL model 4J| it includes the Polyakov loop 
as an effective field and adds quark degrees of freedom. In 
this configuration the EoS describes chiral restoration as 
well as the deconfinement phase transition, while it con- 
tains the correct asymptotic degrees of freedom (quarks 
O hadrons). For details, we refer the reader to Ref. [40| 



B. Thermal contributions to the dimuon excess: 
emission rates 

During the locally equilibrated hydrodynamical stage, 
dimuon emission is calculated locally in space-time, i.e. 
for each cell of the (3-1-1) hydrodynamical grid, accord- 
ing to the expression for the thermal equilibrium rate of 
dilepton emission per four-volume and four-momentum 
from a bath at temperature T and baryon chemical po- 
tential pb ■ The total dimuon yield is then obtained inte- 
grating the emission rate over all fluid cells and all time 
steps of the (3-)-l) grid that are spanned by the system 
during the hydrodynamical evolution until the transition 
criterium is reached. 

In the low invariant mass region of the dilepton spec- 
trum (M < 1 GeV) the largest contribution to the dilep- 
ton excess is due to the p° — )■ emission. As pointed 
out in various works (see e.g. [ll|, El, EE Ell), when 



invoking the vector meson dominance model, the ther- 
mal rate can be related to the spectral properties of the 
vector mesons. Retaining only the contribution of the p 
meson (which is the largest), one arrives at the following 
expression (46i] : 



d-^xd^q 



M2~ 



fB{qo;T) Im DpiM,q;T, ps) , 

(3) 

where a=e'^/(47r)=l/137 denotes the fine structure con- 
stant, M"^ = Qo ~ q'^ the dilepton invariant mass squared 
with energy go and three- momentum q, fsilolT) the 
thermal Bose distribution function, L{M'^) the lepton 
phase space factor. 



L{M^) = 1 



2mf 



1 - 



Amf 
M2" 



(4) 



that quickly approaches one above threshold, and 
lia Dp(M, q;T, pb) the imaginary part of the in- medium 
p meson propagator 



DpiM,q;T,pB) 



1 



Ep{M,q;T,pB) 



(5) 



Here, mp denotes the pole mass and Ep(M, g; T, pb) the 
in medium self-energy. In this application, the self-energy 
contributions taken into account are the following: 

Sp(M, q; T,pb) = S°(M) + £"^(5; T) + J:^^ (q; T, pb) , 

(6) 

where I]'^(Af) is the vacuum self-energy and I]'"^(g;T) 
and {q;T, Pb) denote the contribution to the self- 
energy due to the direct interactions of the p with, respec- 
tively, pious and nucleons of the surrounding heat bath. 
The latter have been calculated according to Ref. [47{ , 
where they were evaluated in terms of empirical scatter- 
ing amplitudes from resonance dominance at low energies 
and Regge-type behaviour at high energy. In principle, 
the matter part of the self-energy should depend on M 
and q (or on g° and q) s eparately. Here, it depends only 
on q because in Ref. [47{ the scattering amplitudes were 
evaluated on the mass shell of the p meson. 

Our analysis of the invariant mass dilepton spectrum is 
restricted to M<1.5 GeV. In the mass region 1<M<1.5 
GeV the sources expected to give a major contribution to 
the dilepton excess are the multi-pion emission and the 
qq annihilation in the QGP. A third background source 
usually present in the intermediate mass region, the cor- 
related decays of D and D mesons, could be disentangled 
from the prompt excess via vertex reconstruction [48j and 
have been recently subtracted from the NA60 data on the 
excess [3|. 

The four-pion contribution is the major manifestation 
of the multi-pion emission for the region of our interest, 
1<M<1.5 GeV, and will be included in this analysis. 
Six-pion contribution as well as other multi-hadron inter- 
actions are not included since they start to play a non- 
negligible role only for masses Af>1.8 GeV. Following 
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Huang's observation [49j that lepton pair production can 
be determined from the electromagnetic spectral function 



extracted in e+e~ 
bution of the 47r 
as 



annihilation, we estimate the contri- 
— s> l^l^ process to the dilepton rate 



4a" 



Sxd^q (27r)2 



-a{e+e- An) . (7) 



The above expression neglects the soft-pion final state 
interaction corrections that induce the mixing between 
the axial-vector and vector currents. We postpone the 
investigation of the effects caused by the parity mixing 
phenomenon to a future work and, here, restrict ourselves 
to mention that the mixing is expected to enhance the 
rate in the invariant mass region 1<M<1.3 GeV (49| . 
Calculations that include the axial vector-vector mixing 
can be found in Refs. 0, 

For the specific purposes of this work, it is sufficient 
to use an empirical approach and estimate the four-pion 
contribution to the dilepton spectra directly from the 



measured cross sections, without attempting any phe- 
nomenological modelling of the underlying processes. We 
consider the multi-pion reactions e^e^ — ?• 7r+7r^7r+7r^ 
and e~^e~ — >■ Tr+Tr^Tr^Tr". As a side remark, we men- 
tion that models based on effective Lagrangians suggest 
that these reactions are dominated by processes involv- 
ing the two-body intermediate state ai(1260)7r [sol. [5l|. 
For the cross section of the e'^e~ 7r+7r~7r"'"7r~ process 
we use the recent precise BaBar data [S^l which cover a 
large range of center of mass energies . The measured 
cross section is in good agreement with the high preci- 
sion data taken at VEPP-2M by SND [HI and CMD- 
2 [H, HH in the energy range 0.7-1.4 GeV, as well as 
with data obtained at DCI by DM2 ^ in 1.4-2.0 GeV 
range. In a similar spirit, for the cross section of the pro- 
cess e'^e~ — 7r"'"7r~7r'^7r^, we take preliminary results of 
the BaBar Collaboration [57] , which agree with SND [s^l 
measurements in the energy range below 1.4 GeV. 



Finally, the qq — >• 7* 
according to [5^ as 



emission is evaluated 



d'^xd'^q 47r^ q 



fB{qo]T)^el In 



+ exp[-(go + ^iq)/T]) {x+ + exp[-/i^/r]) 
{x+ + exp[-(go + fJ'q)/T]) + exp[-/i,/T]) 



(8) 



with x± = exp[— ((70 ± q)/2T] and /ig the quark chemical 
potential. 

Eqs. dS]), ([7]) and ([HJ are valid in the rest frame of 
the fluid. For a moving fluid the Bose/Boltzmann dis- 
tribution function must be substituted with the Jiittner 
function. The modulus of the three momentum can be 
expressed as a function of the Lorentz invariant vari- 
ables M and q/j.u'^ {u'^ denotes the four-velocity of the 
fluid cell) as -y(g)/u^^p"^^M^, so that, e.g., the in- 
medium p meson propagator in Eq. Q can be expressed 
as Dp{M, y/iqf,ui^y - A/2;r, /is). In brief, the substi- 
tution {qQ,q) — > (qf^uf^, \J {q^^ui^Y ~ ^■^^) must be per- 
formed on the r.h.s. of Eqs. ©, © and ([S]). 
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The x-EoS used in this analysis includes an extended 
crossover between the hadronic and the QGP phase dur- 
ing which the two states of matter coexist. The fraction 
of QGP in the medium is estimated as the ratio between 
the energy density stored in quarks and gluons degrees 
of freedom and the total energy density. We denote this 
ratio as A = X{T^fiq). Typical values of A as a function 
of temperature T and quark chemical potential /iq are 
shown in Fig. [T] For a given configuration of (T, ^q) of 
a fluid cell, we weight the hadronic radiation rate with 
the function (1 — A) and the qq annihilation rate with 
the function A. Thus, each cell contributes to the total 
dimuon emission rate according to: 



d'^xd'^q 



^[l-\{T,^iq)](^ 



d^N, 



d^N, 



d-^xd^q d'^xd'^q 



dN - 



d'^xd'^q 



(9) 
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Cascade contribution to the dimuon excess 



In the evolution stages that precede or follow the hy- 
drodynamic phase, dimuon emission from the p-meson is 
calculated as in Ref. [H, [59| employing the time integra- 
tion method (often called also "shining method") that 
has since long been applied in the transport description 
of dilepton emission (see e.g. [59l - l62| ). Note that in the 



pre- and post-hydrodynamical stages the particles are the 
explicit degrees of freedom and all their interactions are 
treated explicitly within the cascade transport approach. 
This allows to dynamically account for final state inter- 
actions during the late stage of the reaction and to follow 
the continuous decoupling of the different particles, since 
sequential freezeout of different particle species occurs 
depending on the microscopic reaction rates. For such 
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FIG. 1. (Color online) Fraction of QGP for various values 
of temperature and quark chemical potential. The dashed 
line indicates where the change of the chiral condensate with 
respect to T and fj,q has a maximum while the solid line shows 
the same for the change of the Polyakov loop. See Ref. [iol ] 
for details. 



a minutely detailed microscopic description of the evo- 
lution of the system, however, there is a price to pay, 
namely that the cascade model solely treats collisions 
and decays on the basis of vacuum cross sections and 
decay rates. Thus, eventual residual in-medium modifi- 
cations of the p meson spectral function in this last stage 
will be neglected when adopting the present model. In 
fact, the consistent treatment of broad spectral struc- 
tures in transport approach is not trivial. Man y w orks 
have been dedicated to this topic, e.g. Refs. |63l - [69| : for 
an overview we refer the reader to the recent critical re- 
view by J. Knoll [73. 

Emission from the stage preceding the hydrodynami- 
cal evolution is typically small, since the geometrical cri- 
terium adopted to start the hydrodynamical evolution 
corresponds to a starting time istart ~ 1-16 fm at top 
SPS energy. Emission from the stage that follows the 
hydrodynamical evolution receives two main contribu- 
tions: when merging the hydrodynamical stage to the 
UrQMD model to perform the final cascade, the hydro- 
dynamic fields are mapped to hadrons according to the 
Cooper-Frye equation. At this point a certain number 
of primary p^'s are created and enter the cascade. If 
soon after the transition the system is decoupled with 
respect to processes involving p mesons, during the cas- 
cade these primary mesons simply decay, no further 
p"'s are generated and the corresponding dilepton yield 
is determined by the abundance of p'^ created at the tran- 
sition times the dilepton branching. If the system is not 
decoupled with respect to processes involving p mesons, 
as it is presumably in reality and in the present model 
(as we will show), p meson (re)generation and absorption 
will occur through processes such as tttt annihilation and 
resonance decays. These processes will delay the decou- 
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FIG. 2. (Color online) Charged particle density as a function 
of the impact parameter. 



pling, increase the emission time and, consequently, the 
dilepton yield. 



III. COMPARISON TO NA60 DATA 

A. Centrality selection 

The NA60 Collaboration has recently presented data 
fully corrected for geometrical acceptance and pair ef- 
ficiencies of the NA60 detector 01. The acceptance- 
corrected data correspond to nearly minimum bias col- 
lisions, selecting events with a charged particle density 
dNch/dri>30. In order to select the appropriate impact 
parameter range in our simulations, we first simulate 
minimum bias collisions and determine the charged par- 
ticle density as a function of the impact parameter. The 
result is shown in Fig. [51 We find that dNch/dr]>30 cor- 
responds to b<9 fm. With this selection, we obtain an av- 
erage charged particle density {dNch/dr])=115, value that 
deviates from the measured one {dN ch / dr})=12{) only by 
4%. 



B. Invariant mass spectra 

In Fig. [3]we show results for the invariant mass spectra 
of the excess dimuons in various slices in the transverse 
momentum of the dilepton pair pT- The theoretical spec- 
tra are normalized to the corresponding average number 
of charged particles in an interval of one unit of rapidity 
around mid-rapidity [71]. In the invariant mass region 
Af<0.5 GeV the spectra are dominated by the thermal 
radiation from the in-medium p meson. The pT scaling 
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FIG. 3. (Color online) Acceptance-corrected invariant mass spectra of the excess dimuons in In-In collisions at ISSA GeV for 
various bins of transverse pair momenta. The individual contributions arise from in-medium modified p mesons (dotted-dashed 
line), 47r annihilation (double dotted-dashed line), quark-antiquark annihilation in the QGP (thin full line) and cascade p 
mesons (dashed line). The sum of the various contributions is depicted by the thick full line. Experimental data from Ref. [3]. 
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FIG. 4. Same as Fig. [S] but with the sudden freezeout approximation. 



exhibited compares nicely to the experimental observa- 
tions. 

In the region around the p meson peak, a net domi- 
nance of the cascade contribution is found. This result 
differs quite substantially from the one of previous theo- 



retical work which adopted the sudden freezeout approx- 
imation. It is an intrinsic feature of the hybrid model, 
due to the presence of a final cascade after the hydro- 
dynamical evolution, leading to a continuous and slow 
decoupling of the particles. In fact, the sum of the ther- 
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mal and cascade contribution leads to an overestimation 
of the NA60 data in the peak region for pr<l GeV. With 
increasing pT, however, the agreement with the experi- 
mental data improves. This aspect deserves further inves- 
tigations and discussions but will be postponed to a sepa- 
rated dedicated section. If we adopt the sudden freezeout 
approximation we obtain the result shown in Fig. 21 One 
observes that the yield in the peak region is reduced in 
comparison to results obtained with a full cascade-like 
modelling of the time scale of the freezeout process and 
the global result is comparable to calculations by other 
groups (see Fig.4 of Ref. fs^). The component denoted 
with freezeout p and depicted by the dashed line cor- 
responds to the dimuon yield from decays of primary p 
mesons produced at the transition. 

In the intermediate mass region, l<Af<1.5 GeV, we 
find that emission from the QGP accounts for about half 
of the total radiation. The remaining half is filled by the 
considered hadronic sources. The An annihilation alone is 
comparable to the QGP emission only for A/>1.4 GcV. 
It is lower in the mass region 1<M<1.2 GeV, where, 
however the sum of the thermal and cascade p is not yet 
completely negligible. With exception of the two lowest 
Pt bins, where the data are slightly underestimated in the 
region 1<M<1.2 GeV, the description of the intermedi- 
ate mass region is reasonable. As already highlighted in 
the very low mass region, the pt dependence of the ther- 
mal radiation in the intermediate mass region is compat- 
ible with experimental observations. This latter point 
is noteworthy because in the three thermal approaches 
which have been used so far to compute dimuon spec- 
tra in comparison to NA60 data the transverse expan- 
sion is driven by two free parameters, namely the expan- 
sion velocity and the chosen radial profile (see Table 1 in 
Ref. [t^ for the specific values /profiles used by the three 
groups). In the present model, on the contrary, the trans- 
verse expansion during the hydrodynamical evolution is 
univocally determined by the pressure of the EoS and by 
the initial transverse profile generated by the preceding 
cascade. This represents an important improvement in 
the field of the dynamics of thermal dileptons. 

Finally, we find that the importance of the thermal 
contribution over the non-thermal one decreases with in- 
creasing Pt , in agreement with previous findings ^9ril2i] . 



C. Transverse mass spectra 

We now reverse the analysis presented in the previ- 
ous section and discuss transverse mass spectra of the 
dimuon excess for four different bins of invariant mass. 
Let us first focus on the thermal emission during the hy- 
dro phase and analyse the shape of the transverse mass 
spectra of two sources, the emission from QGP and from 
a hadronic source such as the in-medium p. Dilepton ra- 
diation from these two contributions is shown Fig.[S] For 
each mass bin the corresponding transverse mass spec- 
trum has been separately scaled by an arbitrary factor to 
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FIG. 5. (Color online) Transverse mass spectra of thermal 
dimuon emitted from in-medium p (left panel) and QGP 
(right panel) in four mass windows. The dotted line depicts 
the emission from 47r annihilation processes in the mass win- 
dow 1<M<1.4 GeV. The spectra are arbitrarily normalized 
to increase readability. 



facilitate the visualization of the differences in the slopes 
of the four curves. We observe that the QGP emission 
(right panel) presents similar slopes for all four invariant 
mass bins. Such an independence from the mass clearly 
suggests early emission from a hot low flow source. On 
the contrary, the radiation from the in-medium p (left 
panel) exhibits a continuous hardening with increasing 
mass. Emission from a source with finite flow typically 
leads to such a mass ordering phenomenon. The same 
considerations explain why the 47r contribution, depicted 
by a dotted line in the left panel of Fig. [SJ is slightly 
harder than the thermal p meson contribution in the last 
mass bin. Let us recall Fig. [3] and focus on the invariant 
mass mass region l<Af<1.4 GeV: we note that the Att 
contribution is more copiously localised at masses higher 
than the in-medium p. 

Comparison of transverse mass spectra resulting from 
the present model and NA60 data is shown in Fig. [HI 
The model fails to describe the peculiar rise at low pT 
(mT<0.2 GeV) observed experimentally. At present, no 
consistent physical interpretation of the rise could be 
achieved [3, Q. With exception of the lowest pr re- 
gion, the model compares reasonably the experimental 
data in the lowest and highest mass bins. In the mass 
bin 0.4<Af <0.6 GeV, the overestimation of the yield ob- 
served around the peak region starts to appear, but no 
strong deviation in the shape of the transverse mass spec- 
tra can be observed. This is not the case in the mass bin 
0.6<Af<0.9 GeV, where the most severe overestimation 
of the yield observed at low pT modify substantially the 
resulting shape of the transverse mass spectra and clear 
deviations with respected to the measured spectra are 
observed. 

Finally, the lowest mass bin is dominated by the radi- 
ation from the in-medium vector meson, whereas in the 
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FIG. 6. (Color online) Acceptance-corrected transverse mass 
spectra of the excess dimuons in four mass windows. Experi- 
mental data from Ref. 



FIG. 7. (Color online) Same as Fig. [6l but with the sudden 
freezeout approximation. 



other mass bins the resulting spectra are a composition 
of various sources, with exception of the high px region 
of the 0.6<Af<0.9 GeV bin, where the cascade contribu- 
tion clearly dominates. Turning the discussion around, 
we can state that the lowest mass bin can be considered 
as barometer of the radiation emitted by the in-medium 
p. The reasonable agreement with the experimental ob- 
servations suggests that the latter does indeed follow a 
dynamical path as expected from hydrodynamical mod- 
els. 

If we adopt the sudden freezeout approximation we 
obtain the result shown in Fig. [T] This is practically 
identical to the previous one in the lowest and highest 
mass bins, whereas major differences can be observed in 
the two intermediate mass bins and are caused by the 
reduction of the non-thermal emission. 

To complete our analysis of the transverse mass spectra 
of the dimuon excess, we follow the procedure adopted by 
the NA60 collaboration and perform a quantitative anal- 



ysis of effective slope parameters. The measured spec- 
tra have been divided into several invariant-mass bins, in 
each of which the data have been fitted with the func- 
tion {1 / mT)dN / drriT oc exp(— mT/Toff ), where the effec- 
tive temperature parameter Tes is the inverse slope of 
the distribution. The fit range was taken as 0.4<pT<l-8 
GcV identically to the NA60 fit range @. We apply this 
procedure to the results of our calculation. The resulting 
exponential fits are shown in Fig. [S] In Fig. [9] the in- 
verse slope parameters extracted from the fit procedure 
are plotted vs. the dimuon mass, as done by the experi- 
mental collaboration. We observe, that the model is able 
to qualitatively describe the rise of T^ff with mass up to 
the pole position of the p followed by a drop in the inter- 
mediate mass region observed experimentally. However, 
at the quantitative level, the extracted values of Teff re- 
produce the experimental ones only in the first and last 
mass bin, whereas underestimate them in the second and 
third mass bin. One should note that a variation in the 
fit range, taking e.g. the fit range chosen in Ref. fl^ . 
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would increase the temperature by 8% in the first three 
mass bins and 3% in the latter, so that the fall down 
would persist. However, such a fit range would not be 
consistent with the experimental procedure. 

At present, we can only speculate on the origin of the 
discrepancies in the values of the effective temperatures 
extracted in the mass region 0.4<M<0.9 GeV. The ex- 
citation function of the mean transverse mass of various 
hadrons has been examined in previous works [13, IH, 113] 
where it was shown that, whereas the hybrid approach 
reasonably reproduces the general behaviour of the mean 
Tot as a function of energy for various particle species, 
slightly higher values than the experimental data were 
observed for pions, protons and negatively charged kaons 
produced in nucleus-nucleus collisions at the top SPS en- 
ergies. Similar considerations can be drawn from the 
comparison of the inverse slope parameters extracted 
for the r], uj and <j) mesons by the NA60 Collaboration 
from the analysis of the cocktail (Fig. [9l open circles) 
with the respective values obtained within the present 
model (Fig. [HI open triangles). For completeness, the 
model result for the effective temperature of negatively 
charged pions 73] has been inserted in Fig. [3] as well. We 
note that the effective temperatures of the three cock- 
tail hadrons are slightly overestimated by the hybrid ap- 
proach. Thus, a more copious radial flow than the one 
resulting from the present model can be hardly concili- 
ated with the overall hadron results at the same beam 
energy. 

The model predicts for the (p an only slightly higher 
value of Tgff than for the u, in qualitative agreement 
with the experimental observations, but on the same time 
returns comparable values for the w and the p mesons, 
where with p here we intend the contribution to the total 
emission in the mass bin 0.6<M<0.9 GeV from freeze- 
out p as resulting when adopting the sudden freeze-out 
approximation in the calculation. On the contrary, the 
inverse slope parameter Tf.f / extracted by the NA60 data 
for the peak region, isolated from the continuum in the 
mass spectrum by a side- window subtraction method and 
interpreted as the freeze-out p [Sj , is about 50 MeV higher 
than that extracted for the w mesons, of similar mass, 
and above the one measured for protons in Pb-Pb colli- 
sions, which are also strongly coupled to pions [tI]. Such 
a behaviour can be hardly obtained in models in which 
decoupling from the hydrodynamical description happens 
simultaneously for all particles species, as the present. 
Nevertheless, as will be shown in the next section, a close 
analysis of the cascade stage that follows the hydrody- 
namical evolution does suggest that after the transition 
the p meson is still coupled to the medium. Within the 
present model, however, essentially no additional radial 
flow is developed in this late stage since the cascade de- 
scription is characterized by softer transverse dynamics 
than the hydrodynamics one, so that the maximum cou- 
pling to the flow is typically reached at the end of the 
hydrodynamical evolution. This may suggest that the 
stage during which the p emerges from the hydrodynam- 



ical flow, with corresponding modifications of its spectral 
function, is too roughly described by the approximation 
that, instantaneously and within a layer, one may transit 
from an in-medium modified p coupled to the fluid to a 
vacuum-like p interacting according to a cascade descrip- 
tion. Presumably, residual in-medium modifications are 
still present during the time-span of decoupling and the 
cascade stage. Escape probability arguments, as well as 
the analysis of the cascade stage, suggest that the ability 
to decouple is momentum dependent since fast particles 
escape earlier than slow ones. To take all these effects 
into account properly is not trivial and beyond of the 
aims of the present work. Generally such models can be 
constructed as e.g. done by Grassi et al. [ll]. Indeed, 
an indirect evidence for the presence of in-medium mod- 
ifications of the low pt cocktail w's has been observed 
by the NA60 collaboration in the form of disappearance 
of the yield in the low tot region with respect to a ref- 
erence exponential fit line @. If one assumes a similar 
behaviour for the p mesons, then a depletion at low rriT of 
the freeze-out contribution to the transverse mass spectra 
can be expected. Whether this can effectively result in 
larger values of the inverse slope parameter of the whole 
excess depends substantially on the relative importance 
of the freeze-out and thermal contribution in the low tot 
region. 

Another effect that may harden the transverse dynam- 
ics is the inclusion of viscosity in the hydrodynamical 
equations, as schematically shown in Refs. [7^,[73]. How- 
ever, the viscous effects induce a hardening of the trans- 
verse spectra not only of dileptons but also of pions (tJ , 
therefore a consistent interpretation of the transverse 
spectra of hadrons and dileptons at top SPS might be 
challenging even for viscous hydrodynamics. 



IV. DECOMPOSITION OF THE CASCADE 
CONTRIBUTION AND DISCUSSION 

In this section we want to have a closer look to the 
various components that enter the contribution of the 
cascade stage to the dimuon spectra. The total cascade 
emission consists of: (i) the emission from the stage that 
precedes the hydrodynamical evolution (denoted as "pre- 
hydro"); (ii) the dimuon emission from the primary p's 
produced at the point of transition between the hydro- 
dynamic and the cascade evolution via the Cooper-Frye 
equation (denoted as "Cooper"); (iii) the emission from 
p mesons (re)generated from interactions occurring dur- 
ing the late cascade stage that follows the hydrodynamics 
evolution (denoted as "regenerated"). 

In Fig. [TO] we show the three contributions for three 
different bins on transverse momentum, a low (left), an 
intermediate (center) and a high one (right). The emis- 
sion from the stage that precedes the hydrodynamical 
evolution (dotted line) is negligible in comparison to the 
others. The contribution to dimuon emission from p 
mesons generated in the cascade stage (dotted-dashed 
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FIG. 8. (Color online) Exponential fit of the transverse mass spectra of the excess dimuons in four mass windows. The fit 
range is restricted to the transverse momentum interval 0.4<pT<1.8 GeV (Note that the fit range in mr — M is different). 



line) is in general comparable to the one from p mesons 
created at the transition (dashed line), with its relative 
importance decreasing for increasing transverse momen- 
tum. This latter pattern is due to the generally more 
moderate transverse dynamics of cascade approaches in 
comparison to hydrodynamical models, as pointed out 
by various studies on transverse mass spectra of different 
hadron species performed within the present framework 
[23I l36l . Isrj and by a previous independent comparative 
analysis of hadronic cascade dyiiamics vs. hydrodynam- 
ics done by Huovinen et al. [78|. A more explicit view 
of the effect is shown in Fig. [TT] where the transverse 
mass spectra of the three components are presented sep- 
arately. We see that the contribution to dimuon emission 
from p mesons that are merged into the cascade at the 
transition point via the Cooper-Frye equation and emit 
in the cascade stage of the evolution (dashed line) has 
flatter spectra than the one from p mesons generated 
in the cascade stage (dotted-dashed line). The emission 
from the stage that precedes the hydrodynamical evolu- 
tion presents as well steeper spectra than the one from 
p mesons generated at the transition point (the latter is 
maximally coupled to the flow). The reason is that this 
contribution is mainly due to early decay of p mesons 



generated by string excitation in first chance nucleon- 
nucleon collisions and therefore reminds closely the steep 
trend exhibited by dimuon transverse spectra in pp col- 
lisions. 

Going back to the comparison with NA60 data shown 
in Fig. [3J we found that the sum of the thermal and the 
cascade contributions lead to an overestimation of the 
peak region. The reason for the disagreement might be 
twofold: 

• As mentioned in the previous section, it can be ar- 
gued that in medium modifications of the spectral 
shape of the p meson are still important during 
the cascade stage and cannot be neglected as typi- 
cally done in cascade models. The authors of Ref. 
[T2I], e.g., took partially into account in- medium 
modifications of the freezeout contribution by em- 
ploying, at freezeout, the vacuum form of the p 
self-energy augmented with a width corresponding 
to the full-width-half-maximum of the in-medium 
spectral function at the freezeout conditions. This 
procedure allows to effectively account for reso- 
nance decays figuring into the p self energy with 
the net effect of depleting the yield in the peak 
region. On the contrary, if resonance decays are 
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and open squares, full circles) are from Ref. 0, The values 
of Tefi extracted by the fit of the calculated transverse mass 
spectra (shown in Fig. (6]) are indicated by open diamonds. 
The full diamonds depict the values of Tofi extracted by the 
fit of the transverse mass spectra obtained using the sudden- 
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composition into peak and continuum of the p-like window 
(see Ref. [sS] for details). The corresponding hybrid model 
results are indicated by open triangles. 



simply added channel by channel in a perturbative 
way, the net effect is the opposite, i.e. to enhance 
the yield in the peak region (this property is quite 
general, see e.g. discussion in Ref. llOl)- It 
might be interesting to investigate these two sce- 
narios more in detail in the future. 

It can be as well argued that the modelling of the 
freezeout time scale via hadronic cascade delivers 
a too slow decoupling of the p meson as a conse- 
quence of a longer persistence of the "pion wind" 
via processes such as e.g. tttt — >■ p ^ tttt and 
ttN N*/A* ^ pN. Concerning this point, 
independent information might be obtained from 
HBT analyses. The analysis of HBT radii of pions 
produced in heavy-ion collisions at the SPS energy 
regime performed within the hybrid approach used 
in this work indeed suggested that at top SPS a 
shorter duration time of the freezeout process than 
the one obtained from the inclusion of a cascade 
stage after the hydrodynamical evolution may be 
favoured by experimental data [79| . An HBT anal- 
ysis of dilepton emission might help to shed some 
light on the decoupling time of the p meson and 
presumably reveal interesting features on the space 
extension of the sources contributing to the various 
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invariant mass regions. 



V. SUMMARY AND CONCLUSIONS 

In this work we employed an integrated Boltz- 
mann-|-hydrodynamics hybrid approach based on the Ul- 
trarelativistic Quantum Molecular Dynamics transport 
model with an intermediate (3-1-1) dimensional hydrody- 
namic stage to analyze the hadronic contribution to the 
dimuon excess observed in In-|-In collisions at ii^iab=158A 
Ge"V. This is the first time that dilepton emission both 
in the thermal and non-thermal regime is studied within 
such a macro-|-micro hybrid approach. We found that 
three regions can be identified in the dilepton invariant 
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spectra. The very low mass region of the spectrum is 
dominated by thermal radiation, the region around the 
p meson peak is dominated by late stage cascade dilep- 
ton emission and the intermediate region receives both 
contributions from hadronic and QGP emission, with the 
QGP accounting for about half of the total emission. The 
invariant mass regions A/<0.5 GeV and Af>l, dominated 
by thermal radiation, show reasonable pt scaling reveal- 
ing that if thermal rates are convoluted in a dynamical 
model a comprehensive interpretation of M and pr path 
observed experimentally can be achieved. The model, 
however, fails in describing the region around the vector 
meson peak for low pT, mostly due to the presence of 
a copious emission during the cascade which follows the 
hydrodynamical evolution. The comparison to experi- 
mental data, seems to disfavour the presence of a long 
lived cascade emission in which the p meson can be ap- 
proximated by its vacuum properties. 

The present calculation represents a first effort towards 
a more consistent treatment of the dynamics of ther- 
mal dileptons from in-medium modified hadrons as well 
as non-hadronic sources, which was so far entrusted to 
more schematic models. In this paper, we discussed the 



dimuon invariant mass region 2m^<M<\.b GeV. The 
NA60 experiment, however, has succeeded in isolating 
the excess also in the region 1.5<M<2.6 GeV [il], and 
a nice matching between the low and the intermediate 
mass region analyses was achieved lis, 801. Extension of 
the present approach in order to enable for the investiga- 
tion of invariant masses up to 2.6 GeV is highly desirable 
and planned for the future. We expect it to bring fur- 
ther insights into the role of partonic contributions to the 
dimuon excess measured in the intermediate mass region. 

ACKNOWLEDGMENTS 

The authors acknowledge useful discussions with 
R. Rapp, H. van Hees and B. Bauchle. G. Moschelli 
and T. Kolleger are thanked for their assistence in the 
use of ROOT. We thank S. Damjanovic for providing 
the experimental data. This work was supported by 
the Hessen Initiative for Excellence (LOEWE) through 
the Helmholtz International Center for FAIR (HIC for 
FAIR). Computational resources have been provided by 
the Center for Scientific Computing. 



[1] 

[2] 

[3] 
[4] 

[5] 

[6] 

[7] 

[8] 

[9] 

[10] 

[11] 
[12] 

[13] 
[14] 
[15] 
[16] 
[17] 
[18] 



R. 



Arnaldi 



et 



al. 



(NA60), 
(NA60), 



[Phys. Rev. Lett. 96, 162302 (2006)| 

R. Arnaldi et al. 

[Phys. Rev. Lett. 100, 022302 (2008)J 

R. Arnaldi et al. (NA60), Eur. Ph ys. J. C61, 711 (20097] 

K. Arnaldi et al. (NA6OJ, 

|Phys. Rev. Lett. 10 2, 222301 (2009) | 

G. A gakishie v et al. (CERES), 
|Phys. Rev. Lett. 75, 1272 (1995)'. 

M. A. Mazzoni (HELIOS/3.), 

|Nucl. Phys. A566, 95c (1994)j 

H. van Hees and R. Rapp, 

|Phys. Rev. Lett. 97, 102301 (2006) | 

K. Dusling, D. Teaney, and L Zahed, 

|Phys. Rev. C75, 024908 (2007)[ 
T. Renk and J. 



Ruppert, 
Zahed, 



|Phys. Rev. C77, 024907 (2008) 
K. Dusling and I. 

i Phys. Rev. C80, 014902 (2009)] 
J. Ruppert, C. Gale, T. Renk, P. Lichard, and J. I. 
Kapusta, |Phys. Rev. Lett. 100, 162301 (2008), 
H. van Hees and R. Rapp, 

Nucl. Phys. A806, 339 (2008) 

P. Huovinen, P. V. Ruuskanen, and J. SoUfrank, 
|Nucl. Phys. A650, 227 (1999)| 

P. Huovinen and M. Prakash, 

|Phys. Lett . B450, 15 (1999) | 

P. Huovinen, M. Belkacem, P. J. Ellis, and J. L Kapusta, 
Phys. Rev. C66, 014903 (2002) 

P. Huovinen, P. V. Ruuskanen, and S. S. Rasanen, 
,Phys. Lett. B535 , 109 (2002) 

M. M.' Aggajwal e t al. (WA98), 
|Phys. Rev. Lett. 85, 3595 (2000)' 

A. Dumitru, S. A. Bass, M. Bleicher, H. Stoecker, and 



[19 
[20; 

[21 

[22; 

[23; 

[24; 

[25; 

[26; 

[27 
[28' 
[29; 

[30 
[31 

[32; 

[33; 

[34 

[35; 

[36 

[37 



W. Greiner, |Phys. Lett. B 460, 411 (1999) . 
S. A. Bass et al.. 



Phys. Rev. C60, 021902 (1999) 



T. Hirano, U. W. Heinz, D. Kharzeev, R. Lacey, and 
Y. Nara, "Phys. Lett. B636, 29 9 (2006)| 
B. Baeuchle and M. Bleicher, 

Phys. Rev. C81, 044904 (2010) 

H. Petersen, J. Steinheimer, G. Burau, M. Bleicher, and 



H. Stocker, Phys. Rev. C78, 044901 (2008) 
E 



Santini, H. Petersen, 
| Phys. Lett. B687, 320 (2010) 
B. Bauchle and 

Ph ys. Rev. C 82, 064901 (2010) 
B. 



and M. Bleicher, 



Bauchle 



and 



M. 



M. 



Bleicher, 



Bleicher, 



Phys. Lett. B695, 489 (2011)', 

E. Santini and M. Bleicher, J. Phys. Conf. Ser. 270 
012040 (2011). 

S. A. Bass et a/.,|Prog. Part. Nucl. Phys. 41, 255 (1998) 



M. Bleicher et al, J. Phys. G25, 1859 (1999)^ 
H. Petersen, M. Bleicher, S. A. Bass, and H. Stocker, 

(2008), arXiv:0805.056 7 [hep-phj] 

J. Steinheimer et aL, |Phys. Rev7^C77, 034901 (2008)| 
H. Petersen and M. Bleicher, 



and J. A. Maruhn, 



'Phys. Rev. C79, 054904 (2009) 
D. H. Rischke, S. Bernard, 
Nucl. Phys. A59 5, 346 (19 95) 

D. H. Rischke, Y. Pursun, and J. A. Maruhn, 
|Nucl. Phys. A595, 383 ("T995j | 

H. Petersen, J. Steinheimer, G. Burau, and M. Bleicher, 
Nucl. Phys. A830, 283c (2009) 

J. Steinheimer, V. Dexheimer, H. Petersen, M. Bleicher, 
S. Schramm, et al, Phys.R ey. C81, 044913 (2010) 
Q.-f. Li, J. Steinheimer, H. Petersen, M. Bleicher, and 
H. Stocker, Phys. Lett. B674, 111 (2009)' 
H. Petersen, J. Steinheimer, M. Bleicher, and H. Stocker, 



14 



J. Phys. G36, 055104 (2009) , 
H. Petersen, J. Steinheimer, G . Burau 
|Eur. Phys. J. C62, 31 (2009) 



and M. Bleicher, 



F. Cooper and G. Frye,|Phys. Re v. DIP, 18 6 (1974) , 
J. Steinheimer, S. Schramm, and H. Stocker, 
J. Phys. G38, 035001 (2011) 

P. Papazoglou, S. Schramm, J. Schaffner- 
Bielich, H. Stoecker, and W. Greiner, 

Phys. Rev. C57, 2576 (1998)^ 

P. Papazoglou et al, Phy s. Rev. C59, 411 (1999) 
V. Dexheimer and S. Schramm, 

[Astrophys. J. 683, 943 (2008)' 
K. Fukushima, Phys. Lett. B591, 277 (2004), 
C. Gale and J. I. Kapusta, Nucl. Phys. B357, 65 (199 1)1. 
R. Rapp and J. Wambach, 



Rapp and 

|Adv. Nuc l. Phys. 25, 1 (2000)| 



V. L. Eletsky, M. Belkacem, P. J. Ellis, and J. I. Kapusta, 
|Phys. Rev.'C64 ,'03 5202 (2001)] 
R. Arnaldi et al. 



Arnaldi 
lEur.Phys.J. C59, 



et 

607 (2009) 



(NA60 Collaboration), 



Z. Huang, [Phys. Lett. B361, 131 (1995 )' 

C. Song, C. M. End C. 

Phys .Rev. D50, 1827 (1994)! 

P. Lichard and J. Juran, Phys. Rev. 
B. Aubert et 

Phys. Rev. D71, 052001 (2005) 



Gale, 



D76, 094030 (2007)J 
(BABAR), 



al 



M. 



N. 



Acha 



et 



|J. Exp. Theor. Phys. 96, 789 (2003) 



R 



R. 



Akhmetshin 
Phys. Lett. B475, 190 (2000)" 
Akhmetshin 
101 (2004) 
et 



et 



R 



R. 



et 



al. 



al. 



|Phys. Lett. 



D 



B595, 

Bisello 



al. 

Tl99T)l 



al., 

(CMD-2), 
(CMD-2), 
(DM2), 



|Nucl. Phys. Proc. Suppl. 21, 111 

V. P. Druzhinin, (2007), arXiv:0710.3455 [hep-exj' 
J. Cleymans, J. Fingberg, and K. Redlich, 
Phys. Rev. D35, 2153 (1987) 

K. Schmidt et al., Phys. Rev. C79, 064908 (200 9) [ 
G.-Q. Li and C. M'. Ko, Nucl. Phys. A582, 731 ( 1995)| 
S. Vogel et al.,^Pi^s. Re v. C78, 044909 (2008)^ 
E. L. Bratkovskaya and W. Gassing, 

Nucl. Phys. A807, 214 (2008)| 

[63] J. KnoU, Prog. Part. Nucl. Phys. 42, 177 (1999)[ 



[64] 
[65] 

[66] 

[67] 

[68] 

[69] 

[70] 

[71] 
[72] 

[73] 

[74] 

[75] 

[76] 
[77] 

[78] 

[79] 

[80] 



S. Leupold, Nucl. Phys. A672, 475 (2000) 
Y. B. Ivanov, J. Knoll, and D. N. Voskresensky, 

S. Juchem, 

D. N. Voskresensky, 



and C. Greiner, 



and C. Greiner, 



|Nucl. Phys. A672, 313 (2000)| 
W! Gassing and 

|Nucl. Phys. A677, 445 (2000) | 
J. Knoll, Y. B. Ivanov, and 
Annals P hys. 293, 126 (2001) 
S. Juchem, W. Gassing, 
|Phys. Rev. D69, 0250 06 (2004) 
S. Juchem, W. Gassing, 
Nucl. Phys. A743, 92 (2004) 

B. Friman, C. Hohne, J. Knoll, S. Leupold, J. Randrup, 
R. Rapp, and P. Senger, eds.. The CBM Physics Book: 
Compressed Baryonic Matter in Laboratory Experiments, 
Lecture Notes in Physics, Vol. 814 (Springer, 2011). 
The rapidity window covered by the NA60 detector is 
about one unit in the forward direction. Here, we make 
use of the symmetry of the system and, in order to in- 
crease statistics, we apply the specular cut |i/|<l around 
mid-rapidity. 

R. Rapp, J. Wambach, and H. van Hees, in Relativistic 
Heavy Ion Physics, Landolt-Bornstein: Numerical Data 
and Functional Relationships in Science and Technol- 
ogy - New Series, edited by R. Stock (Springer, 2010) 
arXiv:0901.3289 [hep-ph] 

For the pions, the fit has been performed in the range 
0<mT — mo<0.7 GeV; for the other hadrons, the same 
fit range as used for the excess, i.e. 0.4<pt<1.8 GeV, has 
been used. 

R. A rnaldi et al. (NA60 GoUaboration), 
[J.Phys G G37, 094030 (2010) 

F. Grassi, Y. Hama, and T. Kodama, 

[Phys.Lett. B355 9 (1995)] 

K. Dusling and S. Lin, |Nucl. Phys. A809, 246 (2008)| 

I\ Song, K. G. Han, and G. M. Ko, 

IPhys. Rev. C83, 024904 (2011)| 

P. Huovinen, P. Kolb, U. W. Heinz, P. Ruuskanen, and 
S. Voloshin, Phys.Lett. B503, 58 (2001) 
G. 



Alt 



et 



al. 



GoUaboration), 



'Phys.Rev. C77, 064908 (2008)', 

S. Damjanovic, R. Shahoyan, and H. J. Specht (NA60), 
GERN Gour. 49N9, 31 (2009). 



